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ABSTRACT 


Recently iTsuii fc Nakaiimal (12014J) and iTsuii. Nakaiima. fc Takedal (1201511 
have developed a method of molecular line spectroscopy of M dwarfs with which 
carbon and oxygen abundances are derived respectively from CO and H 2 O lines 
in the K band. They applied this method to Gliese 229A(G1 229A), the primary 
star of the brown dwarf companion, Gliese 229B(G1 229B). The derived abun¬ 
dances of G1 229A are log Ac = —3.27 ± 0.07 and log Aq = —3.10 ± 0.02, which 
are close to the classical values of the solar abundances of carbon and oxygen. We 
generate model spectra of G1 229B for the metallicity of G1 229A as well as for the 
classical solar metallicity (log Ac = —3.40 and logAo = —3.08). We hnd that 
the differences of the resulting spectra are not so large for the differences of the 
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metallicity of 0.1 dex or so, but we now discuss the spectrum of G1 229B on th e 
basis of the reliable metallicity. From the literature flGizis. Reid fc Hawlevll2002h . 
the lower limit to the age of G1 229A is found to be 0.3 Gyr. From the kinematics 
of G1 229A, we evaluate the upper limit to the age of G1 229A to be 3.0 Gyr. 
The observed and model spectra are compared and goodness of £t is obtained 
in the range of model parameters, 750 < < lOOOK, and 4.5 < \ogg < 5.5. 

Among the candidates that satisfy the age constraint, the best combinations of 
model parameters are (Tgfj, log 5 f) = (800K,4.75) and (850K,5.0), while acceptable 
combinations are (750K,4.75),(850K,4.75) and (900K,5.0). 


Subject headings: stars: abundances — stars: brown dwarfs — stars: low mass 
— stars: fundamental parameters — stars: individual(Gl 229A, G1 229B) 


Introduction 


Since its discovery in 1995, Gliese 229B (G1 229B), a compan i on to an M dwarf, 
has b een regarded as the prototype T dwarf flNakaiima et al.l Il995l: lOppenheimer et ah 
19951) . Although many T dwarfs have be en found in the held by large surveys such as 


SPSS (Strauss et ah 


1999 


York et a. 


2000). 2MASS (Burgasser et al.lll999: Skrutskie et ah 


2006). UKIDSS (Warren et ah 


2007 


Lawrence et al.l 120071) and WISE (IMace et ah 


2013 


Wright et al.ll2010[) . G1 229B remains as one of the best studied T dwarfs to date. Although 
the overall energy distribution of G1 229B was reproduce d by early model spectra fairly well 
(iTsuii et al.lllOOGt iMarlev et al.lllOOGt lAllard et al.l 119961) . there were discrepancies in detail 
between the model and observed spectra. By these models, the effective temperature was 
found to be around 900K, but surface gravity and metallicity were poorly known. According 
to the progres s of ob servations in late 1990’s, attempts were made to improve the models. 
Saumon et ah ( 2000 ) showed that (Te{i,\ogg, [M/H]) = (870K, 4.5,-0.5), ('9 40K. 5.0,-0.2) 
and (1 030K, 5.5,0.0) are all acceptable combinations of the parameters, while Iheggett et ah 
(l2002l) claimed that the best £t combination of the model parameters is (Teg, logg, [M/H]) = 
(lOOOK, 3.5,-0.5) which corresponds to an age of 30 Myr. Despite the accumulation of more 
observational data, the determination of stellar parameters is far from conclusive. In the 
mean time, more and more T dwarfs were discovered, and in 2006, a unifi e d nea r-infrared 
spectral classification scheme for T dwarfs was presented bv lBurgasser et al.l (1200611 . Accord¬ 
ing to this scheme, G1 229B was classihed as a T7 peculiar object, and after this classification, 
G1 229B appears to have been forgotten in the literature. 


The metallicity and age of the primary star G1 229A are expected to be the same as 
those of G1 229B. The metallicity of G1 229A is controversial and the derived values so far 
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are [M/H]=+0.15±0 .15 fjMouldllQTSl) . -0.2±0.4 fISchiavon. Barbuv. fc Singhlll997h . and -0.5 
flLeggett et al.l 120021) . The age indicators are the kin ematics of the youn g disk population 
and stellar activity as a flare star, and based on them iLeggett et ahl (120021) claimed that the 
age of G1 229A could be as young as 30 Myr. 

In this paper, we apply newly determined carbon and oxygen abundances of G1 229A 


(iTsuii fc Nakaiimal l20l4 iTsuii. Nakaiima. fc Takedal 120151) to model the atmosphere of G1 


229B. We also reexamine the age of G1 229A from the points of view of the kinematics and 
stellar activity. We then give the combinations of the good fit parameters of G1 229B, which 
satisfy the age constraint. 


2. Carbon and oxygen abundances of G1 229A 


The details of the determinat ion of carbon and oxygen abundances of M dwarfs including 
G1 229A are described elsewhere (jTsuii fc Nakaiimall20l4 ITsuii. Nakaiima. fc Takedall2015[) 
and a brief summary is given here. 


The A-band echelle spectr um of G1 229A with a re solution of 20,000 was obtained at the 


Subaru Telescope, using IRGS flKobavashi et al 


2000l) with t he na tural guide star adaptive 


optics. Before the spectral line analysis, iTsuji fc Nakajimal (120141) calculated the effective 
temperature and the surface gravity. The effective temperature. Teg = 3710K was estimated 
from the newly derived log Teg — M 3.4 relation, where M 3,4 is the absolute magnitude at 3.4 
pm based on the WISE Wl-band flux and the Hipparcos parallax. The surface gravity was 
estimated to be log q = 4.77 usin g the temperature-radius relation and mass-radius relation 
derived by iBovaiian et ahl (120121) . 


The new analysis of spectral lines is based on the fact that pseudo-continua are seen 
both on the observed and model spectra. The observed spectrum of G1 229A is depressed by 
numerous weak lines of H 2 O, by which a pseudo-continuum is created. On the other hand, 
both the true- and pseudo-continuum are evaluated on the model spectrum. The pseudo¬ 
continuum on the model spectrum is generated owing to the recent H 2 O line database. By 
analyzing EWs of the GO lines affected by H 2 O contamination on the pseudo-continua of 
both observed and model spectra, the difficulty of spectral line analysis of the depressed 
observed continuum can be overcome. 


Almost all the carbon atoms are in the form of GO molecules in the M dwarf atmo¬ 
sphere of G1 229A. This situation is not changed by the changes of physical condition in 
the photosphere and this is the reason that the numerous GO lines are excellent abundance 
indicators of carbon. H 2 O molecules consume a large portion of oxygen left after the GO 
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formation and they are stable in the photosphere of G1 229A. Similarly to the determination 
of carbon abundance from CO lines, oxygen abundance is determined from the analysis of 
H 2 O lines. 


The resultant carbon and oxygen abundances are log Ac = —3.27 ± 0.07 and log Aq = 
—3.10 ± 0.02 respectively. The carbon abundance is higher than the classical high solar 
abundance of log Ac = —3.40, while the oxygen abundance is slightly lower than the classical 
high solar abundance of log An = —3.08, and the C/0 ratio of 0.68±0.12 is slightly higher 
than the solar value of 0.48 flGrevesse et al.lIlQQll) . 


3. Brown Dwarf Models of G1 229B 


We compare the observed spectrum of Cl 229B between 0.8 and 5.1 pm obtained from 
the archive by Leggett (http://staff.gemini.edu/~sleggett/LTdata.html) with a set of model 
spectra based on the carbon and oxygen abundances of Cl 229A. The Cl 229B spectrum 
in t he archive is a comp osite of several spectra with flux calibrations based on photome¬ 
try: ISchultz et al.l fl998[) H ST/STIS spectrum, retained their cali bration, attached at 1.023 
pm to Ceballe et ahl ( 1996 ) JHK spectrum fl u x cali brated as per Leggett et ah ( 1999 ) plus 
far-red spectr um fluxed by Colirnowski et ah ( 19981) [HST 0.81 pm and] 1.04 pm photom- 
etry includes lOppenheimer et ahl (119981) L band spectrum recalibrated 2.98—4.15 pm and 
Noll. Ceballe. &: Marleyi (jl997|l 5 pm spectrum calibrated with UKIRT photometry from 
Colirnowski et ah ( 2004 ). 


To compare with the observed spectrum, we first use the model spectra in a dust free sub¬ 
set based on the classical C & O abundances (Ca-series based on case a abundance) from the 
database of the Un ihed Cloudy Models (UCMs: http://www.mtk.ioa.s.u-tokyo.ac.jp/~ttsuji/export/ucm2) 


dTsuiil I2OO2I . 2OO5I) . Also we newly prepare a small grid of model photospheres (without dust 


clouds) for the metallicity scaled to the carbon abundance of Cl 229A, i.e. all the metals are 
increased except for He, Li and O by —3.27— (—3.40) = -1-0.13 dex against the classical solar 
abundances referred to as case a. The oxygen abundance is kept to be —3.10, the value for 
Cl 229A. We use these abundances to generate model spectra of i? = 600 (about the same 
as the resolution of the observed spectrum). In the brown dwarf as cool as Cl 229B, the 
dust layer lies deep in the atmosphere and there is not significant effect on the photospheric 
emission (jTsuiil 120021) . Therefore the use of dust-free models is justified. 


Before the comparison of the observed and model spectra, we hrst investigate the sig- 
nihcance of the effect of the non-solar metallicity of G1 229A on the Cl 229B’s spectrum by 
comparing the classical solar metallicity model (i.e. of log Ac = —3.40 & logAo = —3.08) 
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with the model with the metallicity of G1 229A for Teff=900K and log 5 '= 5 . 0 . Inspection 
of Fignre [T] reveals that CH 4 bands (at 1 . 6 , 2.3, & 3.5 /rm) are slightly stronger while H 2 O 
bands (1.4, 1.8, & 2.7 /im) are slightly weaker in the model based on the metallicity of G1 
229A, as expected for the increased carbon and decreased oxygen. The differences of the 
spectra, however, are not very large at the resolntion as low as 600. In the discnssion below, 
we always nse the models with the metallicity of G1 229A, bnt the differences of the spectra 
based on the classical solar metallicity and on the metallicity of G1 229A shown in Fignre 
[H generally resnlt in the same model parameters, since the effect of metallicity difference 
of 0.1 dex or so on the low resolntion spectra is rather modest. Thns we confirm that the 
abnndance of G1 229B is not mnch different from the classical solar abnnda nce, and this is 


an ad vantage, since the evolntionary models we nse to interpret the spectra flBnrrows et ah 


200 lh are calcnlated only for the solar metallicity. 


The limitations of the model and observed spectra shonld also be noted before their 
comparison. The FT-band methane absorption in the model spectra is calcnlated by a band 
model and the agreement with the observed spectrnm at this portion of the spectrnm is 
not so good. The iF-band methane absorption in the model spectra is calcnlated by a line 
list. However the qnality of the line list is not so great and the agreement of the model and 
observed spectra at this portion (2.2 — 2.4/im) is not so good either. The model spectra are 
based on eqnilibrinm chemi stry and a featnre related to diseqnilibrinm chemistry, the 4.6 
/rm GO fnndamental band flNolh Geballe. &: Marlevlll997t) is not reprodnced. On the other 
hand, the GO absorption featnre in the observed spectrnm is very noisy and it cannot be 
analyzed simnltaneonsly with other featnres. Signal-to-noise ratios at individnal data points 
of the observed spectrnm are nnknown and the test based on the absolnte valne of is 
not possible. 


Becanse of the limitations mentioned above, we focns onr attention to good portions of 
the observed and model spectra, namely, wavelength segments at (1.0 —1.35/im), (1.45 — 1.6), 
(1.8 — 2.2) and the red wing of L-band methane featnre(3.6 — 4.0/rm) for which a line list 
(R. Freedman, private commnnication) is available for modeling. Since we do not know the 
signal-to-noise ratio of each data point, as proxy for a we define a relative rednced y^, 
y^ by assigning a constant signal-to-noise ratio, G = 1 as 


xl = 




N -I 


C 


Jjiobs) 


ifjiobs) - fj{model)) 


( 1 ) 


where j indicates the jth data point and N is the nnmber of data points. Since C = 
fj{obs)/aj, C/fj{obs) = l/<Jj, where aj is the standard deviation, which is nnknown like C. 
There is no significance in the valne of C. Since we are dealing with only one object G1 
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229B, the usefulness of the relative reduced Xr is justihed. There are 1416 data points in all 
and the htting parameter is the radius of the model brown dwarf. 

A wide range of parameter space is explored in the model htting. We calculate 30 model 
spectra with = 750, 800, 850.900.950 and lOOOK and with log^f = 4.5, 4.75, 5.0, 5.25 and 
5.5. The results of model parameters and goodness of ht, Xr cire given in Tabled! From the 
radius obtained by htting and log 5 ^ of a model brown dwarf, its mass is derived. 


The bolometric luminosity of G1 229B has three determinations. Com bining spec¬ 


troscopic and photometric data from 0.82 to 10 um. iMatthews et a 
6.4 X 1 O“®L 0 . With their own JHKL' photometry, Iheggett et al 

6.6 T 0 -6 X 10 ~^Lf 7 ^. From a recalibration using the HST photometry of iGolimowski et al. 


(1199611 found L = 
'fll999|l obtained L = 


( 1998 1. ISaumon et all ( 2000 1 found L = 6 . 2 ±0.55 x 10 ^Lq. We here adopt t he mid value of 
the bo lometric luminosity L = 6.4 x IO-^Lq and the evolutionary models by iBurrows et al.l 


(1200 ill to estimate the ages of the model brown dwarfs, which are given in the column 5 of 
Table dl 


4. Age of G1 229A 
4.1. Kinematics of G1 229A 


The space velocity of G1 229A is {Uq,Vq,Wq) = (± 12 ,- 11 ,- 12 ) in units of kms“^. 
The solar moti on in the Local Standar d of Rest (LSR) frame of reference is (f/, V, W) = 
(10.0,5.2,7.2) ( Dehnen fc Binnev 1998 1 and the velocity of G1 229A in the LSR frame is 
v{LSR) = 25.5 kms“^, which is signihcantly smaller th a n the LSR velo city of the average 
disk population of n = 43 kms“^ (IGilmore fc Zeiliklll999l) . iLeggettI (jl992[l classified this star 
as a young-disk star, which is reasonable from the value of the LSR velocity. 


A quantitative estimate of the age is possible if we assume that the present LSR velocity 
of G1 229A is produced by a stochastic process in the galactic disk. The total velocity 
dispersion (kms“^) in the LSR frame, as a function of stellar age r (Gyr), is described by 


o'viT)^ = 1000 ± 8.19 X 10'^[exp(r/T) — 1], 


( 2 ) 


where T is 5 Gyr (IGilmore fc ZeilikI 1199911 . 


The fraction of stars with random velocities smaller than v{LSR) for a given velocity 
dispersion a, F{a), is given by the isotropic Boltzmann distribution. 




































( 3 ) 


F{a) 


rv{LSR) 



v'^dv. 


The values of F{a) are given for r and a in Table [H According to the table, the 
probability that G1 229A is actually younger than 3.0 Gyr is 1 — F{a) = 94%. We adopt 3.0 
Gyr as the upper limit to the age of G1 229A. 


4.2. Stellar activity of G1 229A 


While G1 229A is shown to flare, the time-integrated energy in flare light is very much 
lower than for ‘normally active’ flare stars. M dwarfs with high chromospheric activity 
show Balmer lines in emission, while G1 229A shows Balmer lines in absorption indicating 
a low level of chromospheric activity flBvrne. Dovle fc Menzie.sl Il985h . The stellar activity 
indicates that G1 229A is no t extremely young, but not so old. As a quantitative analysis, 
Gizis. Reid fc Hawlevi (l2002l) obtained a lower limit to the age of G1 229A, based on the 
absence of Ha emission and their argument is summarized in the following. 


The analysis of M dwarfs in open clusters by iHawlev. Reid fc Tourtellotl fll999[) suggests 
an approach to an age-activity relation. While the activity levels of stars in a given cluster 
exhibit considerable scatter, there is a well-dehned V — Ic color at which activity becomes 
ubiquitous. All stars redder than this color are dMe (dehned as EW(Ha)> I.OA), while 
the bluer stars are dM without emission. Observations exist for M dwarfs in six clusters: 
IG2602 and IG 2391(30 Myr), NGG 2516 and Pl eiades (125 Myr), Hyades (625 Myr), and 
M67 (4.0 Gyr). I Hawley. Reid fc Tourtellotl (1199911 have used these observations to determine 
the relationship between the “Ha limit” color and the age. 


1/-/c = 2.54+ 1.05 log r. 


( 4 ) 


where r is the age in Gyr. For G1 229A. lGizis. Reid fc Hawlevi (1200211 used Ha absorption 

on th is age limit, Gizis 


and V — I = 2.01 to derive the lower limit to the age, 0.32 Gyr. Base c 


et ah rule out the young age of 30 Myr advocated by Iheggett et ah 


200211 and suggest that 


the synthetic spectra from model atmospheres are not yet adequate for age determinations. 
We adopt 0.3 Gyr as the lower limit to the age of G1 229A, taking into account the error in 
the use oiV — I instead oiV — Iq- 

























5. Interpretation of Models of G1 229B 


From the previous section, we have given the range of the age of G1 229A to be between 
0.3 and 3.0 Gyr. In the column 7 of Tabled], whether a model satishes this age constraint or 
not is given by Y or N. The highest gravity models for \ogg = 5.5 are all ruled out, because 
the model brown dwarfs for < 900K are more massive than the hydrogen burning limit 
and the model brown dwarfs for Tes = 950K and lOOOK are older than 3 Gyr. The relative 
reduced Xr is given in the column 6 and the ranking of the goodness of £t(GOF) is given 
by A{xl < 0.070), B(0.070 < < 0.075), and G(0.075 < Xr) i'^ column 8 . 

There are two models which satisfy the age constraint(Y), and GOF=A, whose physical 
parameters are (Tgff, log ( 7 , m, t) = (800K,4.75,0.028MQ,l.lGyr) and (85OK,5.O,O.O37M0,2.3Gyr). 
These are the best £t models of all, and the spectrum for the former model is plotted in 
Figure [21 

There are three models with (age constraint, GOF) = (Y,B). Their physical parame¬ 
ters are {T^g, logg,m,t) = (750K,4.75,0.038MQ,2.5Gyr),(850K,4.75,0.021MQ,0.70Gyr), and 
(900K,5.0,0.028MQ,l.lGyr). We consider that these models are still good £t to the observa¬ 
tion. As an example, the spectrum for the parameters, (900K,5.0,0.028MQ,l.lGyr) is plotted 
in Figure [3l 

The rest of the seven models which satisfy the age constraint are ranked G. There is no 
(Y,A) or (Y,B) models for T^g =950 and lOOOK. The high temperature models apparently 
£t poorly compared to the models for lower temperatures. As an example of the poor £t 
models, the spectrum for parameters,(lOOOK,5.0,O.O16M0,O.48Gyr), is plotted in Figure ID 

To summarize, we consider that the spectra for (Y,A) models are the best hts and those 
for (Y,B) models are reasonable hts to the observed spectrum. In terms of the effective 
temperature and surface gravity, (Tg^, log 5 f) = (800K,4.75) and (850K,5.0) are the best £t 
parameters, while (750K,4.75),(850K,4.75) and (900K,5.0) are still acceptable model param¬ 
eters. 


6. Comparison with Other Works 
6.1. Metallicity of G1 229A 


A s far as we ar e aware, there have been three determinations of the metallicity of G1 
229A. iMouldl (1197811 obtained high resolution spectra of H and K band using the Fourier 
Transform Spectrometer at KPNO with a resolution R ~ 20,000. He analyzed Al, Ga, and 
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Mg lines and found [M/H]=+0.15±0.15. The effective temperature and the surface gravity 
he adopted are Teg = 3614K and log^f = 4.75 re spectively, which are in go od agreement with 
the estimates Teg = 3710K and log^f = 4.77 by Tsuii fc Nakaiima ( 2014 1. 


Schiavon. Barbuv. fc SinghI fllQQTI ) obtained [Fe/H]=-0.2±0.4 by the analysis of the FeH 
Wing-Ford band at around 1 /xm from a medium resolution spectrum of i? ~ 12, 500. The 
effective ter nperature of their best fit model of = 3300K is significantly lower t han the 
estirnate by Tsuii fc Nakaiima! (2014 1. The two determinations of the metallicity ( Mould 
19781: ISchiavon. Barbuv. fc SinghI 1199711 are both consistent with solar metallicity (within 
the error bars). 


Leggett et al.l (1200211 fit with model atmospheres a spectrum covering from 0.8 to 2.5 
/im with a resolution R ~ 800 and a spectrum covering from 1.12 to 1.22 /xm with a res¬ 
olution R ~ 5,800. Best fit model parameters are (Teg,log 5 f, [M/H]) = (3700K, 4.0, -0.5) 
and (3700K, 3.5, -0.7). Although the effective temperature is in good agreement with our 
estimate, log^f and [M/H] are significantl y lower than our val ues log^f = 4.77 and [C/H] = 
-|-0.13±0.07. The metallicity obtained bv iLeggett et al.l (1200211 is definitely lower than solar 
metallicity. The application of evolutionary mo dels implies that the age of G1 229A is as 
young as 30 Myr. As mentioned in Section 4.2, iGizis. Reid fc HawlevI (1200211 rule ou t this 
yo ung age and suggest that the problem is in the model spectra by lAllard et al.l (1200111 used 
by Leggett et ah ( 2002 1. 


6.2. Metallicity of G1 229B 


Saumon et al.l (120001) do not use any information about G1 229A and their abundance 
analysis is based on the bolometric luminosity L of G1 229B alone. Since L is known, a 
determination of the surface gravity g hxes Teg, the radius, the mass, and the age of G1 
229B, as well as the m etallicity. It is assumed that evolutionary models for solar metallicity 
([Burrows et al.l 1199711 are applicable for non solar metallicity. They ht the observed and 
model spectra by eye. They analyze the “red” spectrum from 0.83 to 1 /xm {R ~ 2250), 
J{R ^ 2400), H{R ~ 2100), and K{R 2800) spectra separately, and derive the H 2 O 
abundances, or [0/H]s. They interpret [0/H] as the metallicity [M/H]. No over-all fitting 
from 0.8 to 5 /xm is given. For a given log^f, the optimal [M/H] was obtained out of a 
set of four [M/H]s for different wavelength bands. For three values of log^f, they generated 
three models (ABG) with optimal parameters: (Teg, log( 7 , [M/H]) = A(870K, 4.5,-0.5), 
B(940K,5.0,-0.3) and G(1030K, 5.5,-0.1). They quote the error in abundance to be ±0.1. 
In terms of the metallicity the model G is the closest to our models, but log 77 is out of the 
range. The compromise is the model B for which the counterpart of our model is (Tgg, log 77, 
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[M/H]) =(900K,5.0,+0.13). 


Leggett et al.l (120021) use the 0.84—4.15 /iin spectr a for G1 229B taken from Leggett et al. 
(1999) and also use the 4.5—5.1 /im spectrum from Noll. Geballe. fc Marlev ( 19971) to ht 


with the model spectra by eye. The best ht model parameters are (Tes,, log^f. [M/H]) 
= (lOOOK, 3 .5, -0.5). Although th e metallicity [M/H]=-0.5 is in agreement with the value 
obtained bv ISaumon et al.l (120001) . the surface gravity log 5 '= 3.5 is extremely small. The 
inspection of the best ht model spectra in comparison with the observed spectrum reveals 
signihcant discrepancies in Z, H, and K bands. These discrepancies ar e much greater than 
thos e seen in our htting of ZJHK bands. The H 2 O line list used by lAllard et al.l (120011) 


i s of [Partridge &: Sc 


(IBarber et al. 


2006 


iwenke ( 1997 ). while that by us is of the BT2-HITEMP2010 database 


Rothman et al.l 120101) . However, the ehect of line lists may be minor at 


low resol ution. Higher resolut ion spectra at J, iL, and K bands (ISaumon et al.ll 2000 l) are also 
htted by iLeggett et al.l (120021) with the best ht model spectra obtained from the analysis at 
the low er resolution. Althou gh the H 2 O line list used bvISaumon et al.l ( 120001 ) was the same as 
that bv lAllard et al.l ( 1200 ll) . the ht bv ISaumon et al.l (120001 ) is si gnihcantly bet t er esp ecially 
at H and K bands. We suspect that the log^f = 3.5 adopted by Leggett et ahl (2002 ) is not 
adequate and consider that the range of log^f = 5.0 ± 0.5 adopted bv ISaumon et ahl (120001) 
is more appropriate. 


7. Disequilibrium Chemistry 

Our analysis is based on equilibrium chemistry and we do not analyze the 4.6 /im 


for T dwarfs ( 

Toll. Geballe. & Marlev 

1997' 

Oppenheimer et a 

. 1998: 

Griffith & Yelle 

1999; 

Saumon et al. 

2000 l 

Visscher Sz Moses 

2011 

Zahnle & Marlev 

2014 

) and we briehy mention 


here the ehect of the G1 229B metallicity on the CO abundance estimation. In general, 
higher gravities, cooler temperatures and lower metallicity favor C H 4 vs. CO. Some workers 


1997: VLsscher & Moses 2011). w 

rile others give CO abundances for low ([M/Hl< —0.5) and 

solar metallicitv (Griffith & Yelle 

1999:ISaumon e 

: al. 2000). Since the carbon abundance we 

adopt is [C/H] = +0.13 ( 

Tsuii & Nakaiima 

2014 

), CO abundances derived for solar metal- 


licity are of our interest. ICrifhth fc Yelle! (119991) estimate logXco > —4 for [M/H]= 0 . 0 , 


while ISaumon et al.l (120001) give logXco = —3.5 for [M/H]=— 0 . 1 . These values are signih¬ 
cantly higher than for example —4.3 < logXco < —3.7 obtained bv lNolh Geballe. &: Marlev 

(I1997I) . 





















































































































































8. Concluding remark 


It has been two decades, since the discovery of G1 229B as the hrst genuine brown dwarf 
in 1995. It had been the only known T dwarf, before a small number of held T dwarfs 
were discovered in 1999. As the prototype T dwarf, G1 229B was studied intensively until 
early this century and efforts were made to determine its physical properties. There are 
two approaches in the determination of physical properties; one is to derive Tes, loggf and 
[M/H] all from the model htting of the G1 229B spectrum itself and the other is to derive Tes 
and logg from the model htting of the G1 229B spectrum, but use [M/H] obtained from the 
spectrum of the primary star, G1 229A. Although the latter approach is favorable because of 
the fewer parameters in the model htting of the G1 229B spectrum, it was deemed difficult 
because G1 229A is an M dwarf. It was not until recently that the metallicity (abundances 
of G and O) was determined reliably. We used these newly determined G and O abundances 
of G1 229A in the model htting of the G1 229B spectrum. In addition to the metallicity, the 
ages of the binary components are safely assumed to be the same. We estimated the range 
of the age of G1 229A based on its activity and kinematics. This age constraint narrowed the 
combinations of physical parameters of acceptable models using the evolutionary models of 
brown dwarfs. It took twenty years to settle the question of physical properties of G1 229B, 
although early guesses assuming the solar metallicity were not so bad. 


We thank T. Geballe and S. Leggett for providing us with the G1 229B spectrum. We 
also thank the anonymous referee for helpful comments which signihcantly improved the 
manuscript. 
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X [nm] 

Fig. 1.— Solar metallicity model (Ca-series based on the classical C & O abundance) spec¬ 
trum (red line) and a model spectrum for the metallicity of G1 229A (green line) of G1 229B. 
Model parameters are T^s = 900K and log^f = 5.0. GH 4 bands (at 1.6, 2.3 & 3.5 /rm) are 
slightly stronger while H 2 O bands (at 1.4, 1.8 & 2.7 fim) are slightly weaker in the model 
based on the metallicity of G1 229A, as expected for the increased carbon and decreased 
oxygen. The differences of the spectra, however, are not very large at the resolution as low 
as 600. 
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X [nm] 

Fig. 2.— Observed spectrum (red line) of G1 229B and a model spectrum for the metallicity of 
G1 229A (green line) . Model parameters are (Teg, log^f, m, t) = (8OOK,4.75,O.O28M0,l-lGyr). 
This is one of the two models with the ranking of goodness of £t GOF=A that satisfy the 
age constraint. There are two gaps in the observed spectrum between 2.59 and 2.99 pm and 
between 4.15 and 4.51 pm. 
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X [nm] 

Fig. 3.— Observed spectrum (red line) of G1 229B and a model spectrum for the metallicity 
of G1 229A (green line) . Model parameters are {Tes,logg,m,t) = (900K, 5.0, O.O 28 M 0 , 
l.lGyr). This is one of the three models with the ranking of goodness of £t GOF=B that 
satisfy the age constraint. These models show still reasonable fit. 
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X [nm] 

Fig. 4.— Observed spectrum (red line) of G1 229B and a model spectrum for the metallicity 
of G1 229A (green line) . Model parameters are {Tes,\ogg,m,t) = (lOOOK, 5.0, O.OIOM©, 
0.48Gyr). This is one of the seven models with the ranking of goodness of £t GOF=G that 
satisfy the age constraint. We regard that the fits are poor for these models. 





























Table 1. Model Parameters and goodness of fit 


Tft 

K 

logs 

r 

Rq 

m 

Mq 

t 

Gyr 

a 

Ar 

0.3 < t < 3.0 

GOF'" 

750 

4.50 

0.136 

0.021 

0.70 

0.0764 

Y 

c 


4.75 

0.137 

0.038 

2.5 

0.0716 

Y 

B 


5.00 

0.137 

0.068 

>10 

0.0738 

N 

B 


5.25 

0.137 

0.122 

— 

0.0781 

N 

G 


5.50 

0.137 

0.216 

— 

0.0903 

N 

G 

800 

4.50 

0.116 

0.016 

0.47 

0.0752 

Y 

C 


4.75 

0.117 

0.028 

1.1 

0.0696 

Y 

A 


5.00 

0.118 

0.051 

4.3 

0.0673 

N 

A 


5.25 

0.118 

0.090 

— 

0.0685 

N 

A 


5.50 

0.118 

0.161 

— 

0.0748 

N 

B 

850 

4.50 

0.099 

0.011 

0.20 

0.0802 

N 

C 


4.75 

0.100 

0.021 

0.70 

0.0730 

Y 

B 


5.00 

0.101 

0.037 

2.3 

0.0682 

Y 

A 


5.25 

0.101 

0.066 

>10 

0.0667 

N 

A 


5.50 

0.102 

0.120 

— 

0.0705 

N 

B 

900 

4.50 

0.086 

0.009 

0.10 

0.0889 

N 

C 


4.75 

0.087 

0.016 

0.47 

0.0796 

Y 

G 


5.00 

0.088 

0.028 

1.1 

0.0732 

Y 

B 


5.25 

0.088 

0.050 

4.3 

0.0692 

N 

A 


5.50 

0.089 

0.091 

— 

0.0691 

N 

A 

950 

4.50 

0.075 

0.006 

0.06 

0.1010 

N 

C 


4.75 

0.076 

0.012 

0.24 

0.0900 

N 

C 


5.00 

0.077 

0.022 

0.75 

0.0817 

Y 

G 


5.25 

0.077 

0.038 

2.5 

0.0764 

Y 

C 


5.50 

0.078 

0.070 

>10 

0.0739 

N 

C 

1000 

4.50 

0.066 

0.005 

0.04 

0.1145 

N 

G 


4.75 

0.067 

0.009 

0.13 

0.1019 

N 

C 


5.00 

0.067 

0.016 

0.48 

0.0935 

Y 

C 


5.25 

0.068 

0.030 

1.5 

0.0846 

Y 

G 


5.50 

0.069 

0.055 

4.9 

0.0803 

N 

G 


^Relative reduced 

'^Ranking of goodness of fit based on relative reduced A(<0.07)5 

B(0.07-0.075), and C(>0.075). 
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Table 2: Age, velocity dispersion and fraction of stars with v < v{LSR) 


T 

Gyr 

a 

kms”^ 


0.5 

21 

0.312 

1.0 

27 

0.172 

1.5 

31 

0.121 

2.0 

35 

0.088 

3.0 

41 

0.057 

average disk 

43 

0.050 


'F{a) is the fraction of stars whose velocities generated by a stochastic process are smaller than v{LSR) The 
probability that a star is younger than the age t is given by 1 — F{a). 



